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1. Introduction 

In spite of the remarkable agreement of the Standard Model (SM) with available data from 
high-energy experiments, it is expected to be only a low-energy manifestation of a more 
complete theory at energy scales beyond a TeV. This new TeV-scale physics is expected to 
ameliorate the problems that beset the SM because of the huge discrepancy between the 
electroweak scale and the Planck (or GUT) scale. The most popular candidate for such 
an extension of the SM has been its supersymmetric generalisation which, in its simplest 
form, is the Minimal Supersymmetric Standard Model (MSSM). The gauge structure of the 
MSSM essentially replicates that of the SM but, in the Yukawa sector, in addition to the 
usual Yukawa couplings of the fermions to the Higgs (responsible for the fermion masses), 
other interactions involving squarks or sleptons are possible. 

The relevant part of the superpotential containing the Yukawa interactions involving 
squarks or sleptons is given in terms of the chiral superfields by 

Wrpv = —XijkLiLjEk + X[ji^LiQjDk -I- —X'Xj^UiDjDk + iiiLiH2 (1-1) 

where L {Q) are the left-handed lepton (quark) superfields while E, D, and U contain the 
corresponding right-handed fields, and i, j, k generation indices. A and X' are lepton-number 
(L—) violating, the X” couplings are baryon-number {B—) violating and the last term is 
a L-violating bilinear coupling. The simultaneous existence of the L— and B —violating 
couplings can induce a catastrophically high rate for proton decay and are usually forbidden 
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in the MSSM by invoking a discrete symmetry called i?-parity where R = (— 
where S is the spin of the particle, so that the SM particles have -R = 1, while their 
superpartners have R = —1. However, i?-conservation is too strong a requirement to 
avoid the unwanted proton decay for it can be effectively forbidden assuming that either 
the L-violating or the H-violating couplings in Eq. 1.1 are present, but not both. Limits 


on the R-violating couplings derived from existing experimental information have been 
summarised in Ref. j^]. 

In the presence of i?-violating couplings, the lightest supersymmetric particle (LSP), 
which is usually the neutralino, is not stable and can decay through R-violating modes 
This is in contrast to the R-conserving MSSM where the LSP is stable and this stability 
is a very desirable feature if the LSP were to be a viable dark matter candidate. In the 
R-violating case, the neutralino cannot be a dark matter candidate unless the R-violating 
couplings are very small so as to ensure that the lifetime of the neutralino is much more 
than the age of the universe. The situation can be saved, however, in theories where 
the gravitino (the spin-3/2 superpartner of the graviton) is the lightest supersymmetric 
particle: a circumstance that can be realised very naturally in theories with gauge-mediated 
supersymmetry breaking Q. The light gravitino is long-lived enough to account for dark 
matter (or, at least, the hot component of dark matter) even in the presence of R-violating 
couplings @. 


2. Physics of Light gravitinos 

Even though gravity is naturally incorporated if supersymmetry is realised as a gauge sym¬ 
metry, the gravitational sector is usually irrelevant for collider phenomenology because of 
the feebleness of gravitational interactions. But if supersymmetry is broken spontaneously, 
the gravitino acquires a mass by absorbing the would-be goldstino and in the high-energy 
limit the gravitino has the same interactions as the goldstino Q. These interactions are 
proportional to l/uig and consequently the interactions of the gravitino can become im¬ 
portant for processes at collider energies in the uIq —> 0 limit. The mass of the gravitino is 
related to Fq, the fundamental scale-squared of supersymmetry breaking, by the following 
relation: 

niQ = —. ( 2 . 1 ) 

^ VSMp 

Mp = 2.4 X 10^® GeV is the reduced Planck mass, and using this value one obtains 

Fo 


niQ = 5.9 X 10 


-5 


reV. 


(500GeV)2^" 

Given a lower bound on the value of Fq one can then deduce a lower bound on the mass of 
the gravitino which, in turn, yields a bound on the interactions of the gravitino with the 
SM particles. 

To make these considerations more concrete, we write down the relevant part of the 
supersymmetric Lagrangian containing the gravitino interactions: 


£ = 


1 ^ 


8Mp 




1 


\/2iVLj 




(2.3) 
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where G is the gravitino field, the gaugino field, the corresponding field strength 
and {4>, ip) the scalar and the fermionic components of the chiral supermultiplets. At the 
level of a effective interaction, the spin-3/2 gravitino field can be well described by its 
spin-1/2 goldstino component when it appears as an external state, i.e. 


Using this limit in Eq j2.3| , allows one to compute the decay widths of the process Xi 
X jZG^ for example. These are: 


where 


r(x° - iG) 

r(x° ^ ZG) 


■*^27 
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(2.5) 


= \Niicos6w + Ni2sm6w\‘^ 
k-iZt = lA^iisin^w - A/2 Cos6'wP 

i^iZL = \NiscosP - Ni4sm(3f. ( 2 . 6 ) 


The Nij are the Xi components in standard notation. The neutralino can also decay into 
a gravitino and a neutral Higgs particle and the corresponding expressions for these are 


r(x° - </G) = 




m^o 

Xi 


96tt Mpm?^ 


1 - 


2 

Xi 


where the Higgsino components are given by 


Kj/jO = I A’issina — Ai4Cosa|^ 
= I Aiscosa -|- Aij4sina|^ 
= |Aj3sin/3-I-Aj4Cos/?p. 


(2.7) 


( 2 . 8 ) 


From the above equations, it is easy to convince oneself that the decay modes into the 
photon dominates over the decays into the Z or the neutral Higgs boson because the decays 
into the latter states are phase-space suppressed. If the neutralino is bino-dominated, then 
the branching ratio into a photon and gravitino is nearly 100%. 


3. The CDF anomaly 

CDF has recently presented results on the production of combinations involving at least 
one photon and one lepton (e or /i) in pp collisions at 7® = 1-8 TeV, using 86.34 pb“^ 
of Tevatron 1994-95 data |^. In general the results were consistent with the standard 
model, however 16 photon-lepton events with large were observed, with 7.6 ± 0.7 are 
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expected. Moreover, 11 of these events involved muons (with 4.2 ± 0.5 expected) and only 
5 electrons (with 3.4 ± 0.3 expected), therefore there is a clear asymmetry, which indicates 
the existence of a lepton flavour-violating process involving muons. 

As we proposed in an earlier paper Q, we suggest that the excess can be simply 
understood in terms of smuon resonance production via an L-violating A' coupling^ which 
decays predominantly into a bino-dominated neutralino and a muon, with the neutralino 
further decaying into a photon and a gravitino. The production and decay has been shown 
in the Feynman diagram in Figure |l|. The merit of this model that we proposed is that it is a 
natural explanation of the characteristics of the CDF anomaly: 1) the flavour-dependence 
is a direct consequence of the i?-violating coupling and, 2) the fact that the excess is 
seen in final states involving photons emerges very neatly in the model because the decay 
Xi ^ 7 G dominates overwhelmingly over other decay modes. We note that both i?-parity 
violation and the existence of a very light gravitino are needed to explain the anomaly, in 
our model. Nonetheless, we emphasise that if one has i2-parity violating supersymmetry, a 
light gravitino is preferable from dark matter considerations, as explained in the previous 
section. 

A light gravitino has also been previously invoked |^] to explain the ee^'y^rp CDF 
event 0, detected in searches for anomalous production of missing transverse energy 

> 12 GeV), in events containing two isolated, central photons. The event was explained 
in terms of the i?-conserving production of a pair of selectrons and the subsequent decay of 
each of these selectron into a 7 G final state. It has been shown |K] that this explanation is 
excluded in the framework of the minimal uni-messenger gauge mediated supersymmetry 
breaking (GMSB) model, because of the anomalously large rates for jets -I -7 -|- ^rp events 
predicted by this model. The problem can be traced back to the mass spectrum of the 
uni-messenger models: in this version of GMSB models the charginos and second-lightest 
neutralinos are light which lead to large jets +^+^rp rates not seen in experiment. However, 
in multi-messenger models of GMSB the charginos and the second-lightest neutralinos are 
heavier and one can have a viable explanation of the CDF event in these models which is 
not in conflict with other existing experimental information. For our purposes, again a light 
neutralino of 100 GeV mass and a reasonably light smuon in the mass range of about 150 
GeV is needed but we require all the other supersymmetric particles to be very massive. In 
the present paper, we perform detailed fits to the experimentally measured distributions of 
the anomalous events in order to determine the masses of the lightest neutralino and the 
smuon with the assumption that all the other masses are heavy enough not to be produced 
at the Tevatron. We do not attempt to place our scenario in the context of some specific 
model of GMSB, but point out that this is indeed possible in the case of multi-messenger 
models of GMSB. A detailed model-dependent study is relegated to a later publication. 


4. Constraints 

If the anomalous events seen by the CDF experiment are to be attributed to the produc¬ 
tion of a smuon resonance involving an /^-violating operator, we can ask what the precise 

^ Smuon resonances at hadron colliders have been previously studied in a different context Q. 
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Figure 1: Feynman diagram of resonant smuon production followed by neutralino decay. 


form of this operator is. To get a substantial cross-section for the production of the smuon 
resonance one needs to couple it to valence quarks in the initial state. This observation 
is then sufficient to specify the i?-violating operator to be L 2 Q 1 D 1 corresponding to the 
coupling X' 2 ii- This operator generates the interactions find and D^dd (and charge con¬ 
jugates), along with other supersymmetrised copies involving squarks. Therefore, if we 
invoke this operator to explain the CDF anomaly we will simultaneously predict effects 
in other channels which will manifest itself through the production of either sneutrinos or 
squarks. In our model, since we take the squarks to be heavy, their effects on experimental 
observables will be negligible. On the other hand, the sneutrinos are necessarily relatively 
light and can be produced resonantly and should lead to observable effects in experimental 
situations. In the present paper, we not only analyse the smuon resonance production in 
the context of the CDF anomaly but also provide predictions for both the smuon and the 
sneutrino channels at Run I and Run II of the Tevatron. The smuons, sneutrinos and the 
lightest neutralino are all light enough to be pair-produced through R-conserving channels. 
We also provide predictions for these cross-sections at Run I and Run II. 

For our analysis, we have essentially four parameters at our disposal; the gravitino 
mass, itLq, the bino mass parameter Mi, the smuon mass parmeter, rufi and the R-violating 
coupling A 211 . The coupling, A^n, is constrained from = F(7r —> eh')({7T fj,u) |T^ to 
be < 0.059 X iqq 'gev Hi' constraint involves a squark mass which is large 

in our model. So the constraint from R-j^ for our purposes is not very relevant. However, 
instead of simultaneously fitting the four parameters using the experimental data, we choose 
to work with fixed values of the A 211 and and perform fits in Mi and m^. While the 
production of the smuon resonance is through the R-violating mode, its decay needs to go 
through the R-conserving channel to a neutralino and muon final state. The R-violating 
decay of the slepton is possible but constrained, in principle, by the Tevatron di-jet data 

which exclude a a.B > 1.3 x 10^ pb at 95% C.L. for a resonance mass of 200 GeV. 
However, in practice this does not provide a restrictive bound upon our scenario as long 
as the R-violating coupling is sufficiently small. We also add that the di-jet bound is not 
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very restrictive because it suffers from a huge QCD background. By restricting A 211 to be 
small, we also avoid the possible /^-violating decays of the Xi f^jj or Xi final 

states. With these considerations in mind we choose A 211 = 0.01. The gravitino mass is 
also fixed at 10“^ eV in our fits. We will discuss the effects of varying the gravitino mass 
and the /^-violating coupling later in this paper. 


5. Defining the model 


The supersymmetric model parameters that are relevant for our discussion, are: Mi, M 2 , 
/X, tan/3 and mo, which determine the chargino, neutralino and sfermion masses at low 
energies. Since no other exotic cascade decays at CDF, are observed, we assume that: 

• Charginos and other superparticles (except the slepton and the lightest neutralino) 
are too heavy to be produced at the current energies. 

• The decays of the lightest neutralino to gauge bosons other than the photon are 
coupling and/or phase-space suppressed. 

These considerations constrain the allowed supersymmetric parameter space, which 
as we are going to show, still has some generality within it. To see this, let us look at 
the formulae that give the chargino and neutralino masses and mixings in terms of the 
fundamental supersymmetric parameters. 

The tree-level neutralino mass matrix in the —iW^, basis are the mass 

eigenstates of the matrix 


Y = 


( Ml 0 —mz sin cos /3 

0 M 2 mz sin 6w cos (5 

—mz sin Oyi/ cos (5 mz cos 9w cos (5 0 

y m^ sin sin/3 —m^ cos sin/3 —/i 


mz sin 9y/ sin /3 ^ 
—mz cos 9w sin f3 

0 / 


(5.1) 


and are defined by y/ = with Wj being the unitary matrix which diagonalises 

Y. The respective mixings in the basis ( 7 ,.^’) instead of {B,W 3 ), are given by iVj^ = 
Nji cos 9w + Nj2sm9w, ^j 2 = — iVji sin^w + Nj 2 cos 9w, = Nj^ and = Nj^^. 

Finally, the respective chargino mass matrix in the basis is 


X = 


( M 2 mu/\/2sin/3 
miv\/2cos/3 /i 


(5.2) 


In our work we will not make use of the GUT inspired relation Mi = | tan^ 9wM2, but 
will instead keep Mi and M 2 generic. As we see from the above formulas, a light neutralino 
can arise either via a light Mi, a light M 2 , or a light /x. In the second and third cases 
however, the chargino is also going to be light enough to be seen in cascade decays, which is 
not the case in CDF data. Moreover, constraining the relative masses of Mi and M 2 roughly 
determines the photino component of the lightest neutralino. Under these conditions tan /3 
is expected to play a relatively moderate role. Table contains the lightest chargino and 
neutralino masses, and the magnitude of the photino-component of the lightest neutralino, 
for different model parameters. In this section only, we constrain ourselves to small and 
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Ml 

= 90 GeV,tan /5 

= 4 

Ml = 

120 GeV,tan /3 = 4 

M 2 





\Nn\ 

”"x? 



Wu\ 

200. 

-600. 

91. 

202. 

202. 

0.88 

121. 

202. 

202. 

0.88 

200. 

-200. 

90. 

167. 

170. 

0.85 

118. 

168. 

170. 

0.81 

200. 

200. 

77. 

144. 

135. 

0.71 

99. 

151. 

135. 

0.56 

200. 

600. 

88. 

191. 

191. 

0.86 

118. 

191. 

191. 

0.85 

400. 

-600. 

91. 

396. 

397. 

0.88 

121. 

396. 

397. 

0.88 

400. 

-200. 

90. 

197. 

200. 

0.85 

118. 

198. 

200. 

0.83 

400. 

200. 

79. 

191. 

181. 

0.80 

104. 

196. 

181. 

0.75 

400. 

600. 

88. 

381. 

380. 

0.87 

118. 

381. 

380. 

0.87 

600. 

-600. 

91. 

561. 

563. 

0.88 

121. 

561. 

563. 

0.88 

600. 

-200. 

90. 

200. 

202. 

0.85 

118. 

202. 

202. 

0.84 

600. 

200. 

79. 

201. 

191. 

0.81 

105. 

204. 

191. 

0.78 

600. 

600. 

88. 

534. 

533. 

0.87 

118. 

534. 

533. 

0.87 




Ml = 

120 GeV,tan/3 = 50 




200. 

-600. 

119. 

197. 

197. 

0.87 





200. 

-200. 

110. 

158. 

153. 

0.96 





200. 

200. 

109. 

156. 

150. 

0.67 





200. 

600. 

119. 

196. 

196. 

0.86 





400. 

-600. 

119. 

389. 

389. 

0.87 





400. 

-200. 

112. 

197. 

191. 

0.79 





400. 

200. 

111. 

197. 

190. 

0.78 





400. 

600. 

119. 

388. 

388. 

0.87 





600. 

-600. 

119. 

547. 

547. 

0.87 





600. 

-200. 

112. 

203. 

197. 

0.81 





600. 

200. 

111. 

204. 

196. 

0.80 





600. 

600. 

119. 

545. 

545. 

0.87 






Table 1: Light weak gaugino masses and photino component of the lightest neutralino |-/Vii| for 
various values of tan/3, Mi, fj,. 


intermediate values of M 2 and fj,, since for larger values our requirements are more easily 
fulfilled. 

We see that demanding a lightest neutralino in the range 90 — 120 GeV, with the 
chargino remaining relatively heavy leads to a neutralino mixing in a photino component 
that is significantly constrained, and lies in the range (0.81-0.88). This also holds for larger 
values of M 2 ,// which are not included in the table. 

We use a single slepton mass parameter mj = ruji^e^ = Neglecting small fermion 

mass terms, the tree-level slepton masses are then (for the first two generations) 

= m? -h (^ - sin 6'^)M| cos(2/3) 

'me,flu = mf - sin cos (2/3), 
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(5.3) 


+ 2 cos(2/3) 


with negligible mixing proportional to the electron and muon masses respectively. 


We use the ISASUSY part of the ISAJET7.58 package |W] to generate the spectrum, 
branching ratios and decays of the sparticles. For an example of parameters, we choose 
(in the notation used by ref. |W]) A 211 = 0.01, m 3/2 = 10 ^ eV, tan/3 = 10, = 0, 

/i together with other flavour diagonal soft supersymmetry breaking parameters are set 
to 2000 GeV. We emphasise that this is a representative point in the supersymmetric 
parameter space and not a special choice. Any superparticles except the first two generation 
sleptons, the lightest neutralino and the gravitino do not appear in this analysis because 
they are too heavy to be produced or to contribute to cascade decays in CDF data. 


6. Fitting kinematical distributions 

We now simulate the signal events for the process in Figure ||. The Standard Model 
background is taken from ref. |^. We use HERWIG6.3 |jl^ including parton showering (but 
not including isolation cuts) to calculate cross-sections for single slepton production. The 
slepton mass parameter mj and the bino mass parameter^ Mi are allowed to vary in order 
to see what range of neutralino and slepton masses are preferred by the experiment. 

We simulate the detector by the following: 

• Photons can be detected if they do not have rapidity 1.0 < |r/| < 1.1, |r/| < 0.05. The 
region 0.77 < ?/ < 1.0, 75° < cj) < 90° is also excluded because it is not instrumented. 
If these constraints are satisfied, we assume 81% detection efficiency for the photons. 

• Muons have a 60% detection efficiency if |r/^| < 0.6 or 45% if 0.6 < r]f^ < 1.1 

Rapidity of a particle is defined as r] = — lntan(0/2), where 6 is the longitudinal angle 
between the particle’s momentum and the beam, (p is the transverse angle between the 
particle’s momentum and the x-axis. We also implement the cuts used in the experimental 
analysis to beat the background down: £''r(/r) > 25 GeV, > 25 GeV and > 25 

GeV. Because we do not perform jet reconstruction, we do not perform isolation cuts. 

CDF gave one-dimensional projections in the events for the following kinematic 

variables 


Et = \IpI +pI, 

mi2 = y/Pl.P2 

m^ = eI-pI-pI 

A(/>i2 = (pi — 4>2 

Ht = + Ft( 7) -l- Et{p) 

^All parameters take their quoted values at the electroweak scale. 
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( 6 . 1 ) 








where px,y are the x and y (i.e. transverse) components of the momentum respectively and 
pi ^2 refer to 4-momenta of particles labeled by 1 and 2. 

Once we have a statistically large sample of signal events simulated, we have a predic¬ 
tion for the number of expected signal plus background events in bin z of a distribution; 
^S+B- average background presented in ref. |^]. The Poisson distributed 

probability density function (PDF) of observing events compared to Ng+B 


PS+B i^obs ) ^s+B ) 




n:. 


( 6 . 2 ) 


In any one distribution, ps+B = ^iP^s+B gives the total PDF for that distribution, assum¬ 
ing each bin is uncorrelated to the others. Unfortunately, we are not able to make this 
assumption between different distributions of variables, because they contain data on the 
same events and should contain some level of correlation. Finally, the log likelihood of 
signal plus background In ps+B is calculated and normalised by subtracting the analogous 
log likelihood for the Standard Model background prediction 


-21nL =-2(lnp5+B - Inp^)- (6.3) 

A negative value then indicates that the data favour the signal plus background hypothesis 
over just background. When performing parameter (xj) estimation, one determines the 
equivalent number of a away from the best-fit point (which has parameters Xb) by 

(Afj)^ = —2 (InL(xi) — InL(xb)) . (6.4) 

The value of InL(xi) calculated is then equivalent to a probability which matches the 
number of Aa that our model fits the data better than the Standard Model in conventional 
Gaussian statistics. Thus, the measure of number of Act here is purely a measure of 
probability in tails of PDFs, not a statement about the Gaussian nature of that PDF. 

We now perform a fit to Mi and keeping all other fundamental parameters (but not 
branching ratios) constant. We must perform the fit one time for each different distribution; 
because the distributions come from the same events, they are all correlated to some extent. 
We cannot therefore assume that the distributions are uncorrelated in order to fit more 
than one distribution at a time, i.e. the correlations must be taken into account. While we 
can generate the correlations in the signal events by Monte-Garlo, we do not own multi¬ 
dimensional data on the kinematic variables in the data. We therefore cannot perform a 
fit to more than one distribution at any time. However, it is possible to see to what extent 
the individual fits to each variable are compatible with each other. The fit corresponds 
to maximising the log likelihood obtained from Eq. |6.3| for one of the distributions. In 
Table § we present the best fit points, indicating that light neutralino masses, in the range 
of 67-111 GeV are to be expected. For small Mi the lightest neutralino mass is determined 
by Ml, while M 2 controls the chargino mass. If M 2 is heavy (as assumed here), no cascade 
decays involving charginos are kinematically favoured. 
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Variable 

Ml (GeV) 

- Ml (GeV) 

-2A ln(L) 

Aa 

Eriyi) 

87 

35 

-10.94 

3.31 

Et{i) 

67 

30 

-9.36 

3.06 


104 

47 

-6.09 

2.47 


82 

23 

-9.94 

3.15 

MT{^Td) 

96 

44 

-6.54 

2.56 


99 

24 

-8.81 

2.97 


84 

28 

-6.30 

2.51 


111 

33 

-5.28 

2.30 


97 

26 

-8.87 

2.98 


99 

25 

-9.01 

3.00 

Ht 

72 

33 

-5.18 

2.28 

ARf,^ 

83 

24 

-6.26 

2.50 


Table 2: Separate best fit points for each kinematic variable. We display the bino mass parameter 
Ml, the mass splitting between Mi and the slepton mass parameter mp We display the difference in 
log likelihood between our model and the Standard Model —2A ln(L) and the corresponding number 
of sigma Atr that the model fits the kinematic distribution better than the Standard Model. 


For the rest of this section, we concen¬ 
trate on the best-fit point for because 

this gives the best likelihood out of all the fits. 
In Table we show the percentage of events 
making it through each of the cuts for this 
best-fit point. The table shows that 11.4% of 
sleptons produced end up as detected 
events in CDF. The cross-section of 0.091 pb 
predicts a total of 7.86 events in the 
channel. This is higher than the observed 
excess because the Et{p) distribution itself 
prefers it. 

The relevant branching ratios of the smuon 
for this point are 


cut 

percentage 

detected /r 

52.1 

> 25 GeV 

41.7 

detected Exix) > 25 GeV 

20.8 

detected Et{h) > 25 GeV 

11.4 

cr 

0.091 pb 


Table 3: Percentage of SUSY events that sat¬ 
isfy cumulative cuts for events at CDF, 

Run I for the best-fit point. Events that pass a 
cut on a given row also pass those cuts on rows 
above. The 7 -in-active-region cut is described 
in the text. The cross-section after all cuts is 
displayed on the last row. 


BRiilL ^ xV) = 0.984, 
BR{jlL —>■ ud) = 0.015, 
BR{fiL ilG) = 0.001, 


(6.5) 


with a lifetime of 1 x 10 sec, whereas for the lightest neutralino we have 

BRix^i ^ G 7 ) = 0.975, 

BRixi Ge-e+) = 0.019, 


( 6 . 6 ) 


with a lifetime of 1 x 10 sec. At such small values of A 211 and m^, R-parity violating 
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decays of the lightest neutralino are negligible. The light sparticle masses are 


m^o = 86.8 GeV, = 130.8 GeV, = 104.2 GeV, rnej^^fiR = 129.7 GeV, (6.7) 

whereas we have set all of the other sparticles except for the gravitino to be heavy (around 
2000 GeV), so that they play no role in our analysis. 
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Figure 2: Energy distributions for the events. We show the distributions in (a) lepton Et, 

(b) photon Et, ( c ) and (d) EIt- The solid red histogram is signal plus background for our 
best-fit point, the blue dashed histogram is the Standard Model background and the black points 
(with y/N error-bars imposed) are the observed number of events. 


For this best-fit parameter point, we show the predicted distributions of lepton Et, 
photon Et and ^t Ide histograms of Fig. ^ and compare them with the excess of 
the data over the Standard Model background. Aa is labeled on each plot and is the 
equivalent number of sigma that this best-fit point fits a particular distribution better 
than the Standard Model. Obviously the largest Aa = 3.31 is for the lepton Et, since the 
fit is performed to this variable. But we also see that at this point, the other distributions 
also fit the data better than the Standard Model: all fit the data better than the Standard 
model to, at least, 2a except for the fits to Et{'^) for which Act = 1.94. The photon Et 
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Figure 3: Mass distributions for the events. We show the distributions in (a) the invariant 

mass of the /i -7 pair, (b) transverse mass of (c) transverse mass of j — and (d) transverse 
mass of the — — ^ subsystem. The solid red histogram is signal plus background for our best-fit 

point, the blue dashed histogram is the Standard Model background and the black points (with 
'/N errorbars imposed) are the observed number of events. 


seems to be steeper in the data than in either the Standard Model or in our model and this 
is a feature at other values of In Fig. ^ we show the mass distributions. The data 

seems to indicate a bump extra to the Standard Model at lower values of as shown in 
Fig. The angular distributions in Fig. ^ show that our best-fit point fits the observed 
excess well in events where the 7 and roughly back-to-back in Fig. |^. While Fig. 

does not seem a particularly better fit than the Standard Model by eye, nearly all of the 
difference in In L comes from the last bin, where the Standard Model predicted hardly any 
events. 

To calculate 95% C.L. regions, we scan over the parameters Am = m^— Mi and Mi, 
calculating (u)^ from eq. |6.4| at each point and for each kinematical distribution. The 
95%C'.L. is then given by (cr)^ = —2/WiiLbf + 5.99, where \i\Lbf is the log likelihood 
at the best-fit point of the kinematic variable being examined. The 95% C.L. regions of 
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Figure 4: Transverse angular distributions for events. We show the distributions in the 

transverse angle between (a) (d) (c) (d) The solid red histogram is 

signal plus background for our best-fit point, the blue dashed histogram is the Standard Model 
background and the black points (with \/N errorbars imposed) are the observed number of events. 


Am and Mi for each separate fitted kinematical distribution are displayed in Fig. ||. The 
horizontal region at the bottom of the plot displays the LEP bound from neutralino pair 
production where the neutralinos decay to photons and 


|. We note |n] that analysis 
of LEP2 data at the highest energies should be able to cover the region up to Mi = 100 
GeV or so. The “overlap” region Am ~ 30 — 40 GeV and Mi Ri 90 — 120 is encouragingly 
within the 95% conhdence-level regions for all distributions except for Ej'i'j) (shown as the 
area inside the white line in Fig. |^), which prefers Mi < 90 GeV, below the LEP bound. 
Ideally, a correlated fit would be performed to all distributions simultaneously. Then, the 
significance of not having such a good fit for Exi'j) in the overlap region could be calculated. 
The most constraining variables are Et{h) and which require {Am, Mi) < (50,200) 
GeV and (40,150) GeV respectively. The Mt{]^ti1) region constrains M^o to be less than 
120 GeV. The Ht variable is not plotted because it does not constrain any of the parameter 
space at the 95% G.L. 
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Figure 5: Separate good fit regions for each kinematical distribution. The horizontal region at the 
bottom of the plot displays the LEP bound from neutralino pair production where the neutralinos 
decay to photons and JT^ . The horizontal axis is Am = The kinematical variables 

are (a) (yellow), A(j)^^ (magenta), (blue), (black), (green), (red), 

Et{j) (white line) and (b) (red), (green), (cyan) and AR^^ (blue). 


We now display predictions for various quantities overlaid upon the 95% C.L. region 
from the Et{r) distribution for different values of Mi and Am. For example, in Fig. ^a, 
it is shown that the expected number of detected signal events (including the cuts 

described above) is 3-7 in the 95% C.L. region. Each parameter point specifies a particular 
sneutrino mass by eq. |5^, and the A 211 coupling will also lead to resonant production of 
sneutrinos. The sneutrinos decay dominantly into neutrino and neutralino, leading to a 
signal. We use identical cuts to that used for the channel, except for the cuts 
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Figure 6: Expected number of (a) ^ 7 ^ 2 - and (b) 7 -^ 2 - signal events in Run I data. The dashed 
blue curves show (labeled) contours of number of expected events, the dotted green line shows the 
lower bound coming from LEP2 and the red curve shows the region of good fit Et{^) in the for 
events. 


involving muons. Fig. P) shows that between 5 and 35 events of this nature are expected 
within the 95% C.L. region. Standard Model backgrounds should be small, with dominant 
physics background coming from 7 Z production, where Z —> vi>. We emphasise that 
measuring this interesting channel would provide an independent check on our model. In 
an R-parity conserving channel, a bound on the gravitino mass of > 2.7 x 10“® eV 


comes from the non-observation of signal ™ b)0 data |1£]. They place a bound of 
roughly 10 signal events for minimum Ets of 25 GeV at the 95% C.L. for a luminosity 
of 13 pb“^. This would correspond to an upper bound of around 66 if we scale up to 86 
pb“^, as used here. The DO bound is therefore not restrictive^. The predicted rate anyway 
depends heavily upon the values of the unfitted parameters (see conclusions). 

In Figure we show the range of relevant masses and branching ratios over the best- 
fit region. As Figure w indicates, 0.95 < BR{fl —> /ux?) < 0.99 in the best-fit region, 
thus other decay modes of resonant smuon production ought to be suppressed. Similarly, 
0.91 < BR{xi Gj) < 0.98 thus the competing (lepton and 2 jets) and e^e~G decay 
modes of the resulting neutralino are also suppressed to unobservable levels at Run I. These 
branching ratios are dependent upon A 211 and ttIq. The range of viable right-handed smuon 


mass is 130 < < 210 GeV, as shown in Figure |^c. The lightest neutralino 


mass IS 


approximately equal to Mi and varies up to 170 GeV in the best-fit region, as shown in 
Figure 0 d. 

The approximate level of other processes can be roughly estimated by calculating 
the expected numbers of pairs of sparticles at Run I. Neutralino production is predicted 
to be at an unobservable level, but the light sleptons have a non-negligible number of 
expected pairs produced at Run I. Using 86.34 pb of luminosity, we display the expected 


^Also note that our cuts are completely different to those in the DO analysis. 
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Figure 7: Relevant masses and branching ratios in the best-ht region. The dashed blue lines show 
(labeled) contours of (a) BR{fl (b) BR{xi G 7 ), (c) (GeV), (d) m^o (GeV), the 

dotted green line shows the lower bound coming from LEP2 and the red curve shows the region of 
good ht ETi^J■) in the for events. 


number of selection, selectron-sneutrino and sneutrino pairs produced at Run I in fig. |8|. 
No experimental cuts at all have been applied to these events, so detected numbers of 
these events might be expected to be a factor of 5-10 less than the numbers that are 
displayed in the figure. We can see from Figure ^ that in the 95% C.L. region that 
hts the Et{r) distribution, there are between about 0.1 and 1 selection pairs expected, 
depending upon the actual value of the parameters Am and Mi. The dominant decays of 
the selections is e —> ex?, again followed by x? iG. Thus the selection pairs provide 
the correct signature to describe the ee'y'y^T event recorded by CDF at run I. When 
detector efficiencies taken into account, the expected number of events, while 

being less than one, are nevertheless still much higher than the expected Standard Model 
background. Selection pair production is predicted to be at the same level as smuon pair 
production, since they have approximately equal masses and they are produced via gauge 
interactions. However, the muon detection efficiency is somewhat lower than for electrons. 
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Figure 8: Number of (a) selectron, (b) selectron-sneutrino and (c) sneutrino pairs produced at run 
I of the Tevatron. The dashed blue lines show the expected number of events, the dotted green line 
shows the lower bound coming from LEP2 and the red curve shows the region of good fit ET{fi) in 
the for events. 


so the expected number of detected events from smuon pair production is smaller. 

It is also possible to produce sneutrinos, and selectron-sneutrino pairs (with final state 
677 .^2’) predicted to be at the 0.1-2 event level before cuts. In Figure we see that 
sneutrino pair production is predicted to be at the level of 0.1 to 2 events. This final state 
will manifest itself as 


7. Predictions for Run II 

At Run II of the Tevatron, assuming 2 fb“^ of luminosity, our model can be ruled out or 
verified by again looking for an excess in the channel, with much higher statistics. 

For example, assuming the same cuts and detector efficiencies as in our Run I analysis, we 
expect 193 signal events at Run II for our best-fit point because the cross-section increases 
to 0.096 pb for detected events. This estimate will be subject to change once the relevant 
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Figure 9: Expected number of (a) (b) 7 ^^ signal events in Run II data. The dashed 

blue curves show (labeled) contours of number of expected events, the dotted green line shows the 
lower bound coming from LEP2 and the red curve shows the region of good fit in the for 

events. 


cuts and detector efficiencies for Run II are known. We also expect an excess of 740 events 
in the channel from resonant sneutrino production at the best-fit point. We calculate 
the number of events in the channel and display them in Figure assuming 

2 fb“^ of collected luminosity (assuming the same efficiencies and cuts as used in Run I). 
We see that at least 70 events in the channel are expected and at least 150 in the 

channel. This will be sufficient to measure parameters much more accurately, or rule 
the model out completely. 

If a signal is seen at Run II, the kinematic distributions will determine the viable 

parameter space more accurately than Run I data. Sparticle pair production may also be 
viable, and can provide constraints upon the parameter space. For this reason, we show the 
expected number of slepton pairs produced in Figure IT. Figure |l^ a shows that between 2 
and 30 selectron (and smuon) pairs are expected, between 3 and 50 sneutrino-selectron and 
sneutrino-smuon pairs (Figure |T^), and between 2 and 40 sneutrino pairs (Figure |l^). 
Once efficiencies and cuts are taken into account, this adds up to at most a handful of 
events in each channel. Nevertheless, this would provide independent confirmation for our 
scenario and Standard Model background rates would be still extremely low. 


8. Conclusions 

We have demonstrated that R-parity violating supersymmetry with a light gravitino can 
explain an anomalously high measured cross-section for the channel. It also explains 

features observed in the kinematic variables of the signal events. We have used this infor¬ 
mation to constrain the slepton and neutralino mass parameters in the model. Whereas 
we could not perform a combined fit to all the different kinematic distributions, if Mi ~ 90 
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Figure 10: Expected number of sparticle pairs produced at run II of the Tevatron for 2 fb“^ 
luminosity. The dashed blue lines show the expected number of events for (a) selectron (b) selectron- 
sneutrino and (c) sneutrino pair production, the dotted green line shows the lower bound coming 
from LEP2 and the red curve shows the region of good fit for events. 


and Am^ ~ 30 GeV, all of the distributions are fit well. Ideally a combined fit to all kine- 
matical variables would be performed for our model and a measure of the fit probability 
calculated (along with that of the Standard Model). Such a fit would require simulation 
of the backgrounds as well as knowledge of the multi-dimensional distributions of variables 
rather than the one-dimensional projections available to us. We have seen qualitatively 
that the ee'j'y^x event observed in Run I is fit much better by our model than by the 
Standard Model. 

We chose representative parameters ttIq = 10“^ eV, A 211 = 0.01 and tan/3 = 10. 
Varying the first two parameters does not change the kinematics of the event, merely the 
branching ratios of the decays. Thus the total number of signal events in the relevant 
channel changes, but the kinematic shapes in the signal events remain the same. A higher 
value of uIq decreases the number of neutralinos decaying to a photon and a gravitino, but 
this can be compensated for by increasing A 211 to increase the production cross-section. 
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However, at some value of A 211 , R-parity violating decays will dominate over the gravitino 
decays of the neutralino. We also note that changing tan (3 has the effect of changing the 


relationship between mj and the slepton masses, as eq. 5.3 shows. Thus, different values 
of tan/3 could potentially prefer different regions of Am. 


Higher values of A 211 can also be 
accommodated by increasing both the 
mass of the sleptons and the mass of 
the neutralino to decrease the hard pro¬ 
duction cross-section. These points are 
illustrated in turn in table Q, where the 
parameters are all varied to have in such 
a way as to produce a number of events 
comparable to the best fit value of 7.8. 

Point 1 illustrates that a higher value 
for m^, which gives lower branching ra¬ 
tios of Xi ^ 'yG, can be compensated 
by an increase in A 211 , raising the res¬ 
onant smuon production cross-section. 

The kinematic hts still look favourable: 
point 1 fits the data better than the Stan¬ 
dard Model to Afj = 1.7 — 3.3 depend¬ 
ing upon the variable. However, it does 
seem difficult to raise further be¬ 
cause the required increase in A 211 then 
gives dominant decay modes to be R- 
parity violating, decreasing the number 
of signal events Points 2 and 3 

illustrate the insensitivity to values of 
tan (3. We note in point 4 that a heavy 
neutralino of 200 GeV also can provide enough events by raising A 211 , but that the Et{^) 
distribution is predicted to be too hard compared with the data. Also, the signal bump in 
goes to higher energies, disagreeing somewhat with the data. The other distributions 
seem to fit the data quite well however. Point 4 also shows that the predicted number of 
events does vary with A 211 and m^, and so our prediction for the numbers of these 
events is parameter dependent. 


Quantity 

1 

2 

3 

4 

■^211 

0.02 

0.01 

0.01 

0.03 

(eV) 

0.01 

10"3 

10“^ 

10"3 

tan j3 

10 

20 

5 

10 

Am 

35 

40 

35 

35 

Ml 

87 

87 

87 

200 


8.0 

6.7 

7.9 

10.7 

% t 7 

42 

34 

34 

14 

Aa^ExitL)) 

3.3 

3.1 

3.3 

3.3 

Aa{ET{y)) 

1.8 

1.9 

1.9 

-1.1 

Aa{^T) 

2.2 

2.4 

2.4 

1.0 

Aa{m^^) 

2.6 

2.3 

2.6 

-1.5 


2.1 

2.5 

2.3 

1.9 

Aa{MT{^Tl)) 

2.8 

3.0 

2.9 

0.6 

Aa{MT{]pTllA) 

1.7 

1.6 

2.0 

1.0 

Aa{A(tiy^.y) 

2.0 

2.2 

2.1 

2.1 

Aa{A<j)fj_^) 

2.7 

2.7 

2.8 

2.6 

Aa{A(ji^^f,) 

2.8 

2.7 

2.8 

3.1 

AaiHr) 

1.8 

2.1 

2.1 

1.4 

Aa{ARfj_.y) 

2.3 

2.3 

2.4 

2.7 


Table 4: Examples of 4 different parameter points, 
showing the parameters, the number of predicted 
events -^t 7 ) number of tr that the 

point fits the data better than the Standard Model. 


Run H will provide a definitive test of our scenario, by again looking in the y^'y^T 
and channels. Other final states with small SM backgrounds are expected at the 
few-event level. For example, observation of yLyy^rp and eyyy^rp would 

provide independent confirmation of our scenario. We note, however, that towards the 
higher values of Mi expected, less than one event in each of the pair production channels 
is likely once detector effects and cuts are taken into account. 
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